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ABSTRACT 

Precise simultaneous measurements of the frequencies of the two kiloHertz 
quasi-periodic oscillations (referred in the literature as upper and lower kHz 
QPOs) cast doubts on the validity of the simple beat-frequency interpretation 
and some of the modifications introduced to explain the results of the varying 
frequency difference. A new model explains the variation of the frequency 
difference suggesting that the upper kHz QPO, namely Uh is an upper hybrid 
frequency of the Keplerian oscillator under the influence of the Coriolis force 
and the lower kHz QPO is the Keplerian frequency u^- Such an oscillator 
has two branches characterized by high frequency z/^ (~ 1 kHz) and by low 
frequency ul (~ 50 Hz). The frequency ul depends strongly on the angle 5 
between the normal to the neutron star disk and ft - the angular velocity of 
the magnetosphere surrounding the neutron star. In the lower part of the 
QPO spectrum (~ 10 Hz), this model identifies the frequency of radial viscous 
oscillations Uy (previously called "extra noise component") and the break 
frequency Ub, which is associated with the diffusive process in the transition 
region (the innermost part of the disk). According to this model, all frequencies 
(namely u^, u^, and u^) have specific dependences on uk- This paper focuses 
on the verification of the predicted relations. For the source 4U 1728-34, the 
best theoretical fit is obtained for 5 8.3 ± 1.0°, which is shghtly larger than 
6 = 5.5 ± 0.5" previously found for Sco X-1. In addition, wc show that the 
theoretically derived power law relation oc u^-^^ is consistent with the recent 
observations of other atoll and Z-sources. 

Subject headings: X-rays: bursts-stars: accretion disks-stars: neutron — stars: 
individual: 4U 1728-34, Scorpius X-1 
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1. Introduction 

KiloHertz quasi-periodic oscillations (QPO) in low mass X-ray neutron star (NS) 
binaries (LMXB's) discovered by Rossi X-Ray Timing Explorer (Strohmayer et al. 1996, 
van der Klis et al. 1996, Zhang et al. 1996) made it possible to evaluate quantitatively the 
previously suggested beat- frequency model (Alpar & Shaham 1985). Most of the nineteen 
sources with QPO have twin peaks in the kHz part of the spectrum (lower frequency 
marked by K and the higher one by h in the spectrum of the source Sco X-1 shown in 
Figure 1). Within the beat-frequency interpretation, the upper kHz QPO represents the 
Keplerian frequency (van der Khs et al. 1996) while the lower kHz QPO is produced by 
the beating of this frequency with the NS spin frequency, Vspim which has been identified 
with the frequency of nearly coherent oscillations observed during bursts, Vburst or half that 
(Strohmayer et al. 1996, Miller et al. 1998). According to Mendez & van der Klis (1999), 
"if this interpretation is correct, JS.v equals Vspin and should remain constant. The NS spin 
cannot change by ~ 20% — 30% on time scales of days to months." For a number of sources 
it has been shown that Ai/ — I'h ~ decreases systematically with vk (van der Klis et 
al. 1997, Psaltis et al. 1998, Markwardt, Strohmayer & Swank, 1999, Mendez et al. 1998, 
Mendez & van der Klis 1999). These recent findings have undermined the beat-frequency 
model. The new (two-oscillator) model (Osherovich and Titarchuk 1999a and Titarchuk 
and Osherovich 1999, hereafter OT99a and T099, respectively) assumes that the lower 
frequency vk is the Keplerian frequency. Following Titarchuk, Lapidus & Muslimov (1998 
hereafter TLM), we assume a viscous transition layer in the NS disk (see Figure 1) which 
is bounded by the solid surface of the NS (the inner boundary) and by the first Keplerian 
orbit (the outer boundary). The adjustment of the Keplerian disk to the sub-Keplerian 
inner boundary creates conditions favorable for the formation of a hot blob at the outer 
boundary of the transition layer (TLM). When thrown out into the magnetosphere, the blob 
becomes a Keplerian oscillator under the influence of the Coriolis force with two modes: a) 
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a radial mode with frequency 

^H = [^i+{^/^n" (1) 

where is the angular velocity of the rotating magnetosphere, and b) a mode perpendicular 
to the disk with frequency 

vl = {^/T^){i^K/i^h) sin 5 (2) 

where 5 is the angle between 17 and the vector normal to the plane of the Keplerian 
oscillations (6 << 1). From the observed Uh and uk, using (1), one finds as a function of 

and verifies the existence of the angle S which satisfies equation (2) for all observed 
and Uh. For Sco X-1, it was found that S — 0.5° (OT99a). The second oscillator in the 
two-oscillator model has the frequency of viscous oscillations for which the dependence 
of Uy on is obtained through the dimensional analysis of the radial transport of the 
angular momentum in the transition layer controlled by a viscous force determined by the 
Reynold's number, which, in turn, is related to the mass accretion rate (TLM, T099). The 
solution for i/^ can be presented in the form 

t^v = f{^K)CN (3) 

where fiyx) is a universal function for any NS and the constant reflects the properties 
of the specific source. The diffusive process in the transition layer is characterized by a 
break frequency (see Figure 1) which is related to by a power law (T099): 

^6 oc (4) 

For the source 4U 1728-34, the coefficient of proportionality is found to be 0.041 (T099). 
The angle 5, the constant Cjv and the profile of Q determine the frequency range of v^, 
and z/fe. Some of the QPOs have been described previously within the framework of 
the two-oscillator model (T099). For the first time, we identify all predicted frequency 
branches and compare the inferred main physical parameters for two stars (section 2, 3 and 
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4). In section 5, we extend the comparison of the theoretical relation (4) with observed 
spectra from other atoll sources and Z-sources. The summary and discussion are presented 
in section 6. 

2. Classification of QPOs in Sco X-1 

The frequencies of the observed QPOs for Sco X-1 are plotted in Figure 2 as a function 
of pk- From the observed Uh (upper hybrid frequency marked by asterisks in Figure 2) and 
uk, the profile of 0,{uk) has been calculated and modeled using the theoretically inferred 
magnetic multipole structure of the differentially rotating magnetosphere (OT99a). The 
resulting Q profile is 

n/27r = Co + CiuT + C2yT + (5) 

where C2 = 2{CiC^)^/'^. The constants Cq = no/(27r) = 345 Hz, Ci = -3.29 • IQ-^ 
Hz-^/^ C2 = 1-017 ■ 10"^ Hz-5/3 and C3 = -7.76 • 10"^° Hz'^ have been obtained by the 
least-squares fit with = 37.6/39 (OT99a). Using equation (5) for from formula (2) 
we have calculated for different 5. The best fit to the data presented by open circles in 
Figure 2 has been found for 5 = 5.5". The second harmonic of the low branch, 2i>l ~ 90 
Hz observed by van der Klis et al. 1997 has been also fitted by our theoretical curve in 
Figure 2. We employ the distance determination in the function Lebesgue's space, L2 (e.g. 
Korn & Korn, 1961, §15.2.2) defining the Relative Lebesgue's measure Erl as a root of 
the ratio of the sum of the square difference (between the theoretical and observational 
values) and the sum of square of the theoretical values. Presumably Erl is a better measure 
of the difference between observations and theory than the criteria in the case of the 
actual large uncertainty in determining the QPOs and break frequencies. This uncertainty 
is likely significantly larger than statistical errors because the low frequency features are 
broad and lie atop continuum noise of unknown shape. If the data set is not uniform 
for the better control one can use the second measure: the rms deviation, Erms which is 



-6- 



defined as a ratio of the rms and the mean value of the frequency range, (cf. Bendat & 
Piersol 1971). Typically Srl < Srms'-, for vl in Fig- 2 Erl = ^rms =4.6%. The frequencies 
vl and 2ul are marked as HBO in the original spectra for Sco X-1 QPOs presented by 
Wijnands & van der Khs 1999 (Figure 3 in their paper). In the same figure, the break 
frequency was marked (with question mark) and also an extra noise component was shown 
(z/y in our interpretation). We have also analyzed the archival data obtained by the RXTE 
observatory in the course of Sco X-1 monitoring on May 25-28, 1996. Standard tasks from 
FTOOLS v4.2 were used in our data analysis according to the RXTE Cook Book. To obtain 
a broad band power density spectra (PDS), the original 1/8192 s data were regrouped in 
32 s segments with 1/4096 s time resolution data and each spectrum was calculated over a 
continuous time interval with a constant offset. The PDS obtained were fit over the whole 
range (1/32 - 2048 Hz) using the King and Lorentzian profiles. The Lorentzian components 
were added to the model when their statistical significance was higher than Aa. In order to 
characterize the break in each PDS, we assume such a "break frequency" i/;,, at which the 
power drops to 50% of the low frequency plateau value (in terms of the King model). Our 
confidence in the identification of v^, and (Figs. 1-2) for Sco X-1 is based on two tests. 
Firstly, the observations for the two branches in the lower part of the spectra satisfy our 
theoretical relation (4), namely, i^b = 0.041i/^-^^. As a second test, we use the verification 
that the extra noise component has a frequency which satisfies the theoretically derived 
dependence of on v^. In formula (3), the universal function /(z^x) obtained numerically 
from a solution for radial viscous oscillations at the outer edge of the transition region, in 
fact, can be approximated by a polynomial 

i/y = Cat ■ {aiVK + a2i'K + o^^^k + «4^'x) (6) 

where = 1.5033 ■ lO^^, aa = 1.2289 ■ IQ-^Rz-^ = 1.176 ■ lO^^Hz^^ and 

04 = 1.2289 • lO-^^Hz-3. The best fit for {erl = 8.4%) presented in Figure 2 by a 

solid line (marked viscous) is found for Cn — 9.76. The line for the related ui, with the 
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observational points (marked by filled circles) closely matching the theory {erl — 13.9%) 
is also shown in Figure 2. Notice that in vicinity of uk ~800 Hz the hne for Uy crosses the 
line for ul suggesting the possibility of interactions of two oscillators. 

3. Classification of the QPOs in the Source 4U 1728-34 

The classification of QPOs in the source 4U 1728-34 is presented in Figure 3. The 
profile of Q/271 has been modelled using observations for (Mendez & van der Klis 1999) 
and formula (5) with constants Co = 382 Hz, Ci = -5.475 • 10"^ Hz-V3^ C2 = 1.5132 • 10"^ 
Hz-5/3 and C3 = -1.0492 • IQ-^Hz'^ for = 10-5/8 (OT99a). A survey of the observed 
QPO frequencies in the source 4U 1728-34 was presented in Ford & van der Klis 1998, 
(hereafter FVK). The 100 Hz Lorentzian frequency (presented in the Table 1 there) we 
identify with 2ul. The best fit is found for 5 = 8.3 ± 1.0° using the extrapolation of ^^{uk) 
into the range for uk below 550 Hz and also observed ul for uk > 800 Hz. The previously 
identified (T099) corresponds to the LF Lorentzian frequency in the above mentioned 
Table 1 of FVK and vi, corresponds to the broken PL break frequency for three days March 
1, February 24 and February 18. The data for February 16, as explained in FVK and T099 
have to be treated separately (the observed 4.83 ± 0.69 Hz oscillation does not belong to 
the Uf, branch). For this source, the constant Cat is found to be equal to 7.1. In comparison 
with the Sco X-1 classification (Figure 2), the Uy and ui, branches for the source 4U 1728-34 
have a much wider range and therefore are verified with greater confidence: Erl — 8.5% for 
and Er]^ = 19.5% for Uf,.. The branches and 2^^ are not observed simultaneously, but 
still match the theoretical curves {erl — 21% for ul). 

4. Comparison Between QPO Spectra of Sco X-1 and Atoll Source 4U 1728-34 

Due to the slightly higher for source 4U 1728-34, the resulting Q profile has a 
higher frequency than the fl profile for Sco X-1. For Sco X-1, VLrnax/'^T^ ~ 348 Hz, while 
for source 4U 1728-34, Q,rnaxl'^T^ ~ 415 Hz. But even more important for the values of vl, 
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is the difference in 5 between the two stars. For the source 4U 1728-34, 5 = 8.3 ± 1.0° 
and the 2i/x, branch in Figure 3 changes from 124 Hz {vk — 353 Hz) to 182 Hz {vk = 788 
Hz). The same 2i/i branch for Sco X-1 stays below 100 Hz for all observed due to the 
smaller 5 = 5.5 ± 0.5° and the smaller Q [see formula (2)]. On the other hand, for Sco X-1 
the constant Cn — 9.76 is greater than Cat = 7.1 for source 4U 1728-34. Therefore, the 
frequency of viscous oscillations and the related Uh are higher for Sco X-1 for the same 
uk, as shown in Figures 2 and 3. For the same Reynold's number, the greater Cn implies a 
higher viscosity (T099). 

5. Relation Between the Break Frequencies Vb and Viscous Frequency Py 

The relation between and i/y has been approximated by a power law (4) (T099). 
Recent observations of Wijnands & van der Khs (1999) make it possible to verify formula 
(4) for a number of atoll sources and Z-sources. In Figure 4, the observed Uf, are plotted vs. 

on a log-log scale. The open circles correspond to Z-sources and asterisks represent data 
for atoll sources. The straight theoretical hne confirms that the power law ui, Ri 0.041i/^-^^ 
is correct for atoll sources. When is multiplied by a factor of 4, data for Z-sources match 
the same straight line (closed circles represent Aub for Z-sources in Figure 4). Thus for 
Z-sources 

Ub « OMul-^\ (7) 

It is conceivable that the factor of 4 difference in the proportionality constant is due to the 
dominant role of different eigenmodes of the diffusive process for Z-sources in comparison 
with atoll sources. If this is the case, then relation (7) corresponds to the first mode and 
formula (4) represents the second mode of break frequency (see e.g. Titarchuk 1994 for 
the details of the diffusion theory). In general, the diffusion equation suggests 1^6(71) oc n^, 
where n is the integer number of the corresponding eigenmode. The overall Erl including 
the Z-sources reduced points (closed circles) is 27%. 
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6. Summary and Discussion 

In this Letter we have identified all predicted frequency branches for two stars. Their 
frequencies are compared favorably with the two-oscillator model. For the source 4U 
1728-34, the angle S is found {6 = 8.3 ± 1.0°). This is the third star for which S has been 
inferred from QPO observations. For Sco X-1, 5 = 5.5 ± 0.5° (OT99a) and for source 4U 
1702-42, 5 = 3.9 ± 0.2° (Osherovich & Titarchuk 1999b). The predicted relation (4) has 
been confirmed for Sco X-1 and source 4U 1728-34 and also for a number of other atoll and 
Z-sources. The difference by a factor of 4 between the value of Uh for two groups of stars 
suggests the participation of two modes in the diffusive process in the viscous transition 
layer in the disk. By comparison Sco X-1 has greater viscosity in the transition layer than 
the source 4U1728-34, assuming the same Reynold's number. We stress the identification 
of the lowest of the twin kHz frequencies as a Keplerian frequency in contrast to the 
previous views. The solutions for Uy and Uf, have been found by integration along the radius 
R. The conversion of the resulting profiles into the frequency domain is done in according 
to the formula for the Keplerian frequency uk (e.g. T099, Eq.l). The same is true for 
the modelling of the Q profile. The resulting striking correspondence between the observed 
spectra of QPO (Figures 2 and 3) and the predictions of the two-oscillator model we view as 
confirmation of the correct identification. An extrapolation of O— profile is not expected to 
be accurate as the data are insufficient to allow a reliable estimate of the NS spin frequency. 
The complexity of the QPO phenomena still leaves unsolved problems such as ~ IHz and 
~ GHz QPO phenomena. 

We are grateful to Eric Ford, Mariano Mendez, Rudy Wijnands and Michiel van der 
Klis for the data and Joe Fainberg, Bob Stone and particularly the referee for the fruitful 
suggestions. 
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Fig. 1. — Top: Schematic picture illustrating the idea of the two oscillator model. The 
blob forming at the outer boundary of the transition layer oscillates in the equatorial plane 
(h-mode) and also in the plane which is perpendicular to the disk (L-mode). The arrows 
with the notations (K, b, V) specify the sites of the proposed origin of the appropriate 
frequency. For the purpose of illustration we have enlarged the transition layer relative to 
the neutron star. Bottom: Leahy-normalized power spectrum of the Z-source Sco X-1 in 
units of F X P{F), {rms / mean)^ (see the similar power spectrum figure in Wijnands & van 
der Klis 1999). Data presented by sohd lines were averaged over 4 continuous time intervals 
after the observation started on May 25, 1996 at 05:04:13 UTC. The total integrated time is 
~ 10 ks. The best-fit model spectrum (which consists of superposition of King profile and 
five Lorentzian ones) is shown by dashed line. 

Fig. 2. — Classification of QPO in the Z-source Sco X-1. Solid lines are theoretical curves. 

Fig. 3. — Classification of QPO in the atoll source 4U 1728-34. Sohd hues are theoretical 
curves. 

Fig. 4. — The break frequency vs the viscous frequency. Theoretical curve is presented by 
solid line. Observational points (Wijnands & van der Klis 1999) are presented by asterisks 
for atoll and miUisecond X-ray pulsar sources (4U 1728-34, 4U 0614+09, 4U 1608-52, 4U 
1735-44, 4U 1812, 4U 1705-44, IE 1724-3045, and SAX J1808.4-3658) and by open circles 
for Z-sources (Cyg X-2, GX 17+2, GX 5-1 GX 340+0, Sco X-1). The sohd circles are for 
the Z-sources where the break frequency is increased by factor 4 (see text). 
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